772 
VOLUME 81 PAPER No. 772 


PROCEEDINGS 


| AMERICAN SOCIETY 


CIVIL ENGINEERS 
AUGUST, 1955 


11 DEAD-LOAD AND LIVE-LOAD MOMENTS 
IN SHELLS OF REVOLUTION BUILT-IN 
INTO CYLINDERS 


by John V. Huddleston, J.M. ASCE, and 
Mario G. Salvadori, M. ASCE 


DIVISION 


{Discussion open until December 1, 1955} 


Copyright 1955 by the American Society oF Civit ENGINEERS 
Printed in the United States of America 


ENGINEERING MECHANICS 


Headquarters of the Society 
33 W. 39th St. 
New York 18, N. Y. 


PRICE $0.50 PER COPY 


— 


THIS PAPER 


--represents an effort by the Society to deliver 
technical data direct from the author to the 
reader with the greatest possible speed. To this 
end, it has had none of the usual editing required 
in more formal publication procedures, 


Readers are invited to submit discussion apply - 
ing to current papers. For this paper the final 
date on which a discussion should reach the 
Manager of Technical Publications appears on 
the front cover. 


Those who are planning papers or discussions 
for “Proceedings” will expedite Division and 
Committee action measurably by first studying 
“Publication Procedure for Technical Papers” 
(Proceedings Paper No. 290). For free copies 
of this Paper—describing style, content, and 
format—address the Manager, Technical Publi- 
cations, ASCE. 


Reprints from this publication may be made on 
condition that the full title of paper, name of 
author, page reference, and date of publication 
by the Society are given. 


The Society isnot responsible for any statement 
made or opinion expressed in its publications. 


This paper was published at 1745 S. State Street, 
Ann Arbor, Mich., by the American Society of 
Civil Engineers. Editorial and General Offices 
are at 33 West Thirty-ninth Street, New York 18, 
N.Y. 


ae 
| 
at 
{ 


DEAD-LOAD AND LIVE-LOAD MOMENTS IN SHELLS OF 
REVOLUTION BUILT-IN INTO CYLINDERS 


John V. Huddleston,! J.M. ASCE, and 
Mario G. Salvadori,2 M. ASCE 


INTRODUCTION 


In the investigation of thin shells it is commonly found that only membrane 
stresses are developed in large portions of the shell, but that bending mo- 
ments and transverse shears are necessary in some parts of the shell (e.g., 
near boundaries) to eliminate impossible discontinuities predicted by the 
membrane theory. 

The purpose of this investigation is to facilitate the approximate determi- 
nation of bending stresses at the edge of certain types of shells of revolution 
(usually called “domes”). The shells considered have a vertical axis of rota- 
tion, are loaded symmetrically with respect to that axis, and are built-in into 
circular cylinders with the same axis and radius as that of the shell at the 
point of juncture. This investigation deals with the effect of dead and live 
load, and particularly with the evaluation of bending moments and transverse 
shears caused near the junction of shell and cylinder by such load. The 
standard procedure to accomplish this consists in determining first the dis- 
placements and rotations at the edge of the shell due to membrane stresses. 
To restore continuity, the edge moments and shears must produce equal and 
opposite displacements and rotations, which can be computed with good ap- 
proximation by substituting for the portion of shell near the edge a portion of 
the cylinder tangent to the shell, provided the thickness of shell and cylinder 
be small as compared with their common radius (Geckeler’s assumption). (1)3 


Edge Displacements in Membranes of Revolution Due to Dead Load 


Consider a shell of revolution formed by the rotation of a meridian curve 
y(x) about the y-axis (Fig. 1). Let ry and rg be the principal radii of curva- 
ture at any point (x,y) of the shell, r, being the radius of curvature of the 
surface in the meridian plane (i.e., the radius of curvature of the meridian 
itself) and rg the radius of curvature of the surface in a plane perpendicular 
to the meridian. Call § the angle between the normal to the shell and the 
shell axis, and assume that: 

1) the shell has a vertical tangent at the lower edge, i.e., 


0 


dy 


1. Instructor in Civ. Eng., Columbia Univ., New York, N. Y. 
2. Prof. of Civ. Eng., Columbia Univ., New Yori, N. Y. 
3. Numbers in parentheses refer to the bibliography at the end of the paper. 
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dx 


7 at 


from which it follows that 
dr 


2) the derivative of ry with respect to ~ is zero at the lower edge, i.e., 


adr 
1 ot 


3) The thickness of the shell is essentially constant in the neighborhood 


of the edge 
4) Poisson’s ratio is equal to zero. 
Let the symmetrical dead load acting upon the shell at any height y be p 


per unit area of shell, and express p as a function of p: 


P= p, (4) 


where pp is the value of p at 6 = 5 . (The value of p is obtained at any point 


by multiplying the thickness of the shell at that point by the unit weight of the 


material.) 
It has been shown(2) that, under the conditions specified above, the radial 


displacement wo at the edge 9 = iY (positive outward) and the rotation Wp, at 


the same edge (positive outward, i.e., positive if counter-clockwise at 9 = > 
are given by: 


2 


Pe 
° Eh. r 


radius of the equator (i.e., rg at =F or x at y = 0) 
modulus of elasticity of the material 


hg = thickness of shell in the neighborhood of equator 

To = radius of curvature of meridian at the equator (i.e., ry, at O= >). 
The quantity S (the non-dimensional static moment of the load with respect 

to the shell axis) is defined as follows: 


4. This condition is satisfied by all meridional sections symmetrical about 
the equator (x-axis). 
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aS 
(2) 
pa 
where: 


® 
S = r, sin (7) 


and its evaluation for various meridional curves and thickness variations is 
discussed in a later section of this paper. 


Edge Displacements in Shell of Revolution Due to Moments and Shears 


If a uniformly distributed bending moment Mo per unit length of circum- 
ference and a uniformly distributed shear Qo per unit length of circumference 
are applied to the edge of the shell as shown in Fig. 1, the edge will undergo a 
radial displacement wg (positive inward) and a rotation W g (positive inward). 


To find wg and W gs, for shells having a small value of the ratio hg/a (< +” 


it can be assumed without appreciable error that the displacements due to 
Mo and Q, will be felt only in the vicinity of the boundary and that this part 
of the shell can be analyzed as a portion of a long cylindrical shell of radius 
a.(1) On the basis of these assumptions, wg and WY, are given by: 3 


(8) 


Edge Displacements in Cylinder Due to Moments and Shears 


If equal and opposite edge moments and shears are applied to the cylinder 
as shown in Fig. 1, its edge will undergo a radial displacement we (positive 
outward) and a rotation W¢ (positive inward) given by: 


1 
= 
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1 
Ps 8 
1 
(2/9 mM -Q) (9) 
8 8 
where: 
fs Je, Jah, (10) 
Eh 
8 
(11) 
| 


1 
Y. (26M, + (13) 
2 D 
Calling he the thickness of the cylinder in the neighborhood of the edge, B- 


and De are defined by changing the subscript from “s” to “c” in Eqs. (10 and 
(11), respectively. 


Continuity Between Shell and Cylinder 
Two conditions of continuity arise from the fact that there can be no rela- 
tive displacement or rotation between shell and cylinder at the edge. With the 
symbols of Fig. 2 they are: 


Ww 14 
*% *%, (14) 


Substitution of Eqs. (5), (6), (8), (9), (12) and (13) into (14) and (15) leads to 
two equations in the two unknowns Moy and Q>. Solution of these equations per- 
mits writing M, and Q, in the following from: 


a 


h 5 
° 
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\ where: 
1 h, 
f,(r) | g(r) (18) 
1 h, 
q,* f,(r) g(r) (19) 
and: 
(1 - r@) 
(1 +r°) + 2r“(1 +r) 
| at 


3 
a+r’) 


2 
(1 + (1 + r) 


2 
(1 + 
f,(r) 


(1 + r2)* + (1 +r) 


(3 2”) 


(1 + + 2r°(1 +1) 


h (24) 
h,/ 
Values of mo and qo as functions of r for various values of the parameter 


h 
ifs are given in the graphs of Figs. 3 and 4, respectively. For intermedi- 


ate values of i > a linear interpolation gives the correct result since mp 


and qo are linear functions of this parameter. 
Maximum Moment in Shell or Cylinder 


Assuming again that the part of the shell affected by the edge moments and 
shears behaves as a long cylinder, the bending moment M, per unit length of 
circumference of shell at a distance s, from the edge along the meridian may 
be approximated by: (4) 


= - Ps (By) 


where the functions § and { are defined and tabulated by Timoshenko.(5) The 
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4 (21) 


moment M, has a first maximum at a distance §, defined by the equation: 


tan (7 = (26) 
a, -2 


The maximum value of M,(s,) is given by: 


where: 
m= Zz (B,3, (28) 


Similarly, the bending moment M, per unit length of circumference of cylin- 
der at a distance s,, from the edge is given by: 


M(s.) +a (28) 
Ss" 6 
Pe ec 
and has a first maximum at a distance 5, defined by: 


q, 


tan = (30) 
a, +2 yr 
The maximum value of M,(s,) is given by: 


where 


a, = Pet) + (32) 


To enable the determination of the point of greatest internal moment and the 
value of the moment coefficient (mg or m,) for both shell and cylinder from 
a single graph, let the quantities 8 , 8, and m without subscripts apply to 
either the shell or the cylinder. Then the point of maximum internal moment 
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is given by: 


(Pi) 


q 
= 0.760 = for the shell (34) 
° 


4, 
for the cylinder (35) 


© 


| 


yr m 
The value of the coefficient m for maximum moment is now given by: 
m= - AZ (36 
The ratio u = m/m, thus becomes: 


Fig. 5 contains graphs of 8 s vs. A (Eq. 33), and of U vs. A (Eq. 37). 

Fig. 5a represents positive values of A from 0.1 to 100 (in which range “ 

is always negative) and Fig. 5b negative values of A from 0.1 to 100 (for 
which # is always positive). Once values of m, and q, have been determined 
from Figs. 3 and 4, two values of A (one for the shell and one for the cylin- 
der) can be found from Eqs. (34) and (35). Then Fig. 5 can be used for either 
the shell or the cylinder to determine the location and magnitude of maximum 
internal moment. In the range where y is positive, it is always greater than 
unity and the design moment is of the same sign as M, but greater than M,. 
In the range where yw is negative, it may be greater or less than unity and 
gives a design moment of opposite sign to M,. The second peak moment in 
any case (which may be greater than M, but is always small as compared 
with the first peak moment) is located by increasing 8 & by an amount 7 

and is evaluated by multiplying the first peak moment by -e~ ”, which is 

- 0.0432. 

The approximate values of Mz, obtained by means of Eqs. (27) will be 
acceptable for practical purposes provided the shell does not deviate sub- 
stantially from the tangent cylinder. This deviation may be easily checked 
by noticing that: 


‘ 
where: 
and 
6 
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and that the angle @ subtended in the shell by 5, is given approximately by: 


h 
= 0.76 /*s (38) 


It is seen that the Fer deviation i) depends not only on Jh,/a, but also on 
Since < 7/4 for negative values of and in most practical 
cases A is negative for the shell, the angular deviation may be taken to be: 


< 0.60 (39) 


a 


For the case of spherical shells deviations of the order of 0.10 have been 
shown to introduce errors of less than 1% in the moments.(1) Deviations of 
the order of 0.2 are acceptable for practical purposes. (9) 


Case of Small Edge Moments 


When the value of m, obtained from Fig. 3 is very small, a value of A may 
result large enough to be out of range in Fig. 5. It can be shown, however, 
that for A greater than + 100 or less than - 100 the contribution of the first 
term of the right-hand member of Eq. (28) or of Eq. (32) will always be less 
than 2.1 percent of the second term, and that the first term can therefore be 
neglected in comparison with the second, giving the following approximate 
equations for |>100: 


Z for |A\>100 (40) 


c 


These two equations can be further simplified by noticing that for values of 
A larger than 100 in absolute value, 8 § is very nearly equal to 7/4. Using 
this value in Eqs. (40) and (41), the following equations are obtained for mg 
and m,: 


m= - 0.245 q, (42) 


q, 


0.245 — (43) 
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Live Load Moments 


When the load on the shell is a live load defined by: 


P= p, (44) 


where p is the load per unit area of horizontal projection of the shell and pp 
is the load intensity at the lower edge or at the shell top, the edge displace- 
ment and rotations are given by:(2) 


a 


(46) 


In Eq. (45) S is again the nondimensional static moment of the load with re- 
spect to the shell axis: 


t(x) = ax (47) 


Noticing that in the Eqs. (18) and (19) the terms involving | = derive from 
the rotation y , due to the load, and that W, is zero for the live load of Eq. 


(44), the graphs of Figs. 3 and 4 corresponding to i = = 0 may be used 


directly to obtain m, and q, for the case of the live load. The evaluation of 
the maximum moments is then obtained as was done for the case of the dead 
load. 


Evaluation of S and S;, 


The evaluation of S as defined by Eq. (7) can be accomplished by numerical 
integration in most cases and by exact integration in some instances. Results 
will be presented here for ellipsoidal shells generated by half ellipses of 
various ratios of the minor semi-axis b to major semi-axis a (See Fig. 6), 
including spherical shells, for the cases of constant thickness and of thick- 
ness varying linearly with 9. 

For constant thickness, g(f) = 1 and Eq. (7) becomes, upon proper substitu- 
tion for r; and 


sin 4 
+ cos 9) 


This expression can be integrated exactly and gives: 


5. The bending stresses due to live load depend on the single parameter r and 
were also evaluated directly in Ref. (7). 
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| 
b,2 


1 1+k 


k= /1-(2)° 


Values of S for various values of b/a are given in Table I. 


TABLE I 
VALUES OF S FOR ELLIPSOIDAL SHELIS 
OF CONSTANT THICKNESS 


a 


0.2 
0.3 
0.4 
0.5 
0.6 
0.7 
0.8 
0.9 
1.0 


For the case of linearly varying thickness, g(9) = 2 ~, and Eq. (7) becomes: 


2 sing 4 
= = (5) —, £08 (51) 


[sin” + (=) cos” 
This expression can be integrated by numerical methods, and values obtained 
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where: 
0.587 
0.588 
0.637 
0.690 
0.807 
i 0.870 
0.934 
1.000 


by means of Simpson’s one-third rule and h4-extrapolation 8) are given in 
Table II: 


TABLE II 
VALUES OF S FOR ELLIPSOIBAL SHELLS 


g(f) = = 


bd 


0.2 0.114 
0.3 0.171 
0.4 0.228 
0.5 0.295 
0.6 0,361 
0.7 0.428 
0.8 0.497 
0.9 0.566 
1.0 0.637 


For the case of a thickness varying linearly from some value h, at the 
crown to another (greater) value h, at the edge, i.e., for 


h h 
° 


the quantity S can be computed by: 
h a 1 h " 2 


where S, is taken from the table for constant thickness and S2 from the table 
for linearly varying thickness. The value of S computed by Eq. (53) will al- 
ways be less than unity if hy is less than h,. 

In most practical cases, S will be greater than 0.1, and the graphs of 
Figs. 3 and 4 have been extended to include all values of as from zero tol , 


assuming that S will not be smaller than 0.1. = 
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For any other function g(9), when the meridian is elliptical, Eq. (7) can be 
integrated by expressing it in the form: 


s = (2 sin? 


If g(9) can be evaluated pointwise, Eq. (54) can be integrated numerically. 

For more complicated meridian shapes and thickness variations, including 
those specified only by an architect’s drawing, values of S can be determined 
by numerical or graphical methods, using Eq. (7) as a definition and express- 
ing (if necessary) the quantities under the integral sign in terms of other an- 
gles and distances. 

Figure 7 gives the values of S; as defined by Eq. (47) for some of the most 
commonly encountered live load distributions. 


Example 


As an example of the use of the equations and graphs of this paper, consid- 
er an elliptical shell with a = 40 ft. and b = 20 ft. and with a thickness varying 
linearly from 2 inches at the crown to 4 in. = 1/3 ft. at the edge, the shell be- 
ing built-in < ys a cylinder having a thickness of 4 inches. Since for elliptical 
shells rg = b2/a, in this example r, = 10 ft. and the ratio h/r, takes the value 


1/3 _ 
10 ~ 30’ 
load moments, With b/a = 0.5 Tables I and II give 8) = 0.690 and So = 0.295 
and by Eq. (53): 


S = “ (0.690) + (1 - £) (0.295) = 0.492 


hence the present theory will give good approximation of the dead 


With hg = h, = 1/3 ft., the parameter iM oF — takes the value 0.185. From 


Figs. 3 and 4, with r = “a =1, we pe. mM, = - 0.023 and q, = 0.190, and 


from Eqs. (34) and (35): 


= 0.760 


= - 0,760 6.3 for the cylinder 
71(-0,023) 


= - 6.3 for the shell 


Figures 5(b) and 5(a) give then, respectively: 


(3,3, = 0.60 ; = 2.9; 
2.9(-0.023) = 0.067 
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= 0.96 ; = -1.4; 


n= (-1.4)(-0.023) = + 0.032 


The angular deviation in the shell, by Eq. (39), equals: 


£0.60 x #0 0.219 rod 22.5° 


and is small enough to guaramve an acceptable error in the shell sooty 
Assuming the shell material to be concrete weighing 150 lbs. /tt.3 
dead-load intensity at the edge becomes: 


P, = 150 x 0.333 = 50 ibs./ft.° 


2 
and with ro = o = 10 ft.: 


To 0-492 


= 1,310 ft.lbs./ft. 


P, an, = 50x x 


Hence by Eqs. (16), (27), and (31): 


=-0,023 x 1,310 = - 30.2 ft.1bs./ft. 


~0.067 x 1,310 = - 88.0 ft.1bs./ft. 


+ 0,032 x 1,310 = + 42.0 ft.1bs./ft. 


= 0.362 1/ft. 


we finally obtain: 


- 0.60 
= 0.302 = 1.65 ft. 


- 
ate = 2.65 ft. 


With: 
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Thus the shell must be reinforced with outside steel at the edge to resist a 
moment of 30.2 ft. lbs./ft. and outside steel 1.65 ft. from the edge to resist 
a moment of 88.0 ft. lbs./ft., while the cylinder must have the same outside 
steel as the shell at the edge and inside steel 2.65 ft. from the edge to resist 
a moment of 42.0 ft.Ihs./ft. 
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Fig. 6 


o1/e="s 2i/i=" (a) 
on (D/x-1)°d=d 


02/i="S °q(!) | (Pp) 
K 


0 
(0/x-1) d=d 


< 
a 
| 
~ 
a 
° Q 
Qa 
| 
° 
a 
o 
° 
= 
¢ 
Qa 
> 
: 
& 
7172-21 
| 


AMERICAN SOCIETY OF CIVIL ENGINEERS 


OFFICERS FOR 1955 


PRESIDENT 
WILLIAM ROY GLIDDEN 


VICE-PRESIDENTS 


Term expires October, 1955: Term expires October, 1956: 
ENOCH R. NEEDLES FRANK L. WEAVER 
MASON G. LOCKWOOD LOUIS R. HOWSON 


DIRECTORS 


Term expires October, 1955: Term expires October, 1956: Term expires October, 1957: 
CHARLES B. MOLINEAUX WILLIAM S. LaLONDE, JR. JEWELL M. GARRELTS 
MERCEL J. SHELTON OLIVER W. HARTWELL FREDERICK H. PAULSON 
THOMAS C. SHEDD GEORGE S. RICHARDSON 
LLOYD D. KNAPP SAMUEL B. MORRIS DON M. CORBETT 
GLENN W. HOLCOMB ERNEST W. CARLTON GRAHAM P. WILLOUGHBY 
FRANCIS M. DAWSON RAYMOND F. DAWSON LAWRENCE A. ELSENER 


PAST-PRESIDENTS 
Members of the Board 


WALTER L. HUBER DANIEL V. TERRELL 


EXECUTIVE SECRETARY TREASURER 
WILLIAM H. WISELY CHARLES E. TROUT 


ASSISTANT SECRETARY ASSISTANT TREASURER 
E. L. CHANDLER CARLTON S. PROCTOR 


PROCEEDINGS OF THE SOCIETY 


HAROLD T. LARSEN 
Manager of Technical Publications 


DEFOREST A. MATTESON, JR. PAUL A. PARISI 
Editor of Technical Publications Assoc. Editor of Technical Publications 


COMMITTEE ON PUBLICATIONS 
SAMUEL B. MORRIS, Chairman 
JEWELL M. GARRELTS, Vice-Chairman 
GLENN W. HOLCOMB OLIVER W. HARTWELL 
ERNEST W. CARLTON DON M. CORBETT 


